Octet baryon electromagnetic form factors in 
nuclear medium 



G. Ramalho 1 , K. Tsushima 2 , and A. W. Thomas 2,3 

1 CFTP, Institute) Superior Tecnico, Universidade Tecnica de Lisboa, 
Av. Rovisco Pais, 1049-001 Lisboa, Portugal 

2 CSSM and 3 CoEPP, School of Chemistry and Physics, University of Adelaide, 
Adelaide SA 5005, Australia 



Abstract. We study the octet baryon electromagnetic form factors in nuclear 
matter using the covariant spectator quark model extended to the nuclear matter 
regime. The parameters of the model in vacuum are fixed by the study of the 
octet baryon electromagnetic form factors. In nuclear matter the changes in 
hadron properties are calculated by including the relevant hadron masses and 
the modification of the pion-baryon coupling constants calculated in the quark- 
meson coupling model. In nuclear matter the magnetic form factors of the octet 
baryons are enhanced in the low Q 2 region, while the electric form factors show 
a more rapid variation with Q 2 . The results are compared with the modification 
of the bound proton electromagnetic form factors observed at Jefferson Lab. In 
addition, the corresponding changes for the bound neutron are predicted. 



PACS numbers: 13.40.Gp, 21.65.-f, 14.20.Jn, 12.39.Ki 
1. Introduction 

Whether or not hadrons change their properties in a nuclear medium, has been one 
of the long-standing problems in nuclear physics [TJ [2]. QCD is established as the 
theory of strong interactions and quarks and gluons are the degrees of freedom in 
the QCD Lagrangian. It seems natural that in the strong mean fields which pervade 
nuclear matter the motion of the quarks and gluons inside hadrons should be modified. 
Such changes are what is meant by the nuclear modification of hadron properties and 
studying such effect is clearly central to the understanding of dense matter within 
QCD. 

Recently, strong evidence concerning the modification of nucleon properties in 
a nuclear medium has been reported from the proton electromagnetic form factors 
measured in polarized (e, e'p) scattering on 16 [3] and 4 He [U El [7] at MAMI 
and Jefferson Lab. These experiments measured the double ratio of proton-recoil 
polarization transfer coefficients in the quasi-elastic scattering off nuclei, and the 
results were normalized with respective to the double ratio of hydrogen. The results 
from 4 He strongly suggest the modification of the bound proton electromagnetic 
form factors. Furthermore, the study of neutron properties in the nuclear medium 
in [8j, predicts an enhancement of the same double ratio for the neutron, contrary 
to the suppression observed for the proton. The corresponding experiment, to 
measure the polarization transfer, is planned in the future [5]. Theoretically, there 
are several studies concerning the electromagnetic form factors of the nucleon in 
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the nuclear medium [5J [TO] EJ H21 E3] • They are based either on quark degrees of 
freedom [SJ [TUJ US HI], or meson and nucleon degrees of freedom Q2]. However, it is 
very difficult to separate and identify the observed effects in terms of these degrees of 
freedom. In particular, to distinguish a possible change in the nucleon properties in a 
nuclear medium from those of the conventional many-body effects, such as final state 
interactions and meson exchange current, are very difficult, and seems only possible 
in a model dependent way, where experimental measurement involves all such effects 
including the one-body current modification [2l [14j [T5] . Thus, the interpretation of the 
modification observed is still under discussion and has not been settled yet. In these 
circumstances it is helpful to examine the modification of the electromagnetic form 
factors of the nucleon and other baryons in the nuclear medium within alternative 
approaches. 

In this article we study the medium modification of the the octet baryon 
electromagnetic form factors in nuclear matter focusing on the valence quark structure 
of a baryon. Thus, we do not include final state interactions nor meson exchange 
current, where the latter may possibly be important for the magnetic form factors of 
baryons in a nucleus and nuclear matter [13] HU US] [TT1 [HI [19] . For this purpose, we 
use the covariant spectator quark model [20] [21] [22J [23] [24] , which has its basis in the 
covariant spectator theory [25] . The model has been successfully applied to study the 
electromagnetic properties of the octet [3SM2Z] anddecuplet [351 I2H1 1301 1311 132] baryons. 
The model was also very successful in the studies of 7*7V -> A(1232) [33 ] l34l |35] [36 ] . 
7*iV -> 7V(1440) [37 , and 7*iV -> 7V(1535) [38] reactions. For the meson cloud 
effects, we include the pion cloud effects, which are expected to be dominant, in a 
phenomenological manner based on the method applied in (24] [26] [27] . 

In [57] the model was extended to the lattice regime to utilize the lattice QCD 
simulation data, where the electromagnetic form factors of the octet baryons were able 
to be calculated at lattice hadron masses corresponding to large pion mass values used 
in the lattice QCD simulations. Similarly, it is also possible to extend the model to 
the in-medium regime, once we are able to calculate the in-medium modified masses of 
the baryons and mesons appearing in the model. This is the working hypothesis used 
to extend the model, and to calculate the in-medium modifications of the octet baryon 
electromagnetic form factors. For the in-medium masses of the baryons and mesons, 
we use the quark- meson coupling (QMC) model (39] [40], which has been successfully 
applied to study the properties of nuclei [41] 02], hypernuclei [43] [44] [45] [46], and 
hadron properties in a nuclear medium |47j . based on the relativistic valence quark 
structure of hadrons in a nuclear medium. 

Another point to note concerning this study is that the parametrization of the pion 
cloud contributions in vacuum has been improved over what was used in the past |27j 
based on lattice QCD simulation and chiral perturbation theory. In particular, care is 
taken for the neutron charge form factor and the charge radius, where the pion cloud 
contributions are very important. This is also true for the electric charge neutral 
particles. 

This paper is organized as follows. We start by defining the electromagnetic form 
factors in medium in section [2] In section [3] we explain the covariant spectator quark 
model, and describe the electromagnetic currents for octet baryons in the model. The 
extension of the model to the in- medium regime is discussed in section [5] Results are 
presented in section [5] and discussions and summary are given in section [7] 
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2. Electromagnetic form factors in vacuum and in the nuclear matter 

A spin 1/2 baryon B, a member of the baryon octet, has a Dirac structure and 
therefore its electromagnetic structure can be expressed in terms of two independent 
form factors, namely, the electric Geb and magnetic Gmb form factors in vacuum 
(mass Mb) and those in the nuclear medium (mass M B ) G* EB and G* MB , respectively, 
and they are defined below. 



2.1. In vacuum 

Let us consider an octet baryon B, with mass Mb in vacuum. When the initial 
(momentum P_) and the final (momentum P + ) states are on-shell, the electromagnetic 
current (coupling of the baryon with a photon) can be represented as 

r B = p 1B (Q 2 )r + ^b(q 2 )|^, (2.1) 

where q = P + — P_, and Fib and F 2 b are respectively the Dirac and Pauli form 
factors which are the functions of Q 2 = —q 2 . 

Suppressed in equation (12.1[) are the initial (ub) and final state (ub) Dirac spinors, 
functions of P± and the spin projections. For simplicity we represent the current in 
units e = V^bra, with a ~ 1/137, the electromagnetic fine structure constant. 

At Q 2 = they are normalized as 

Pis(0) = e B , F 2B {0) = k b , (2.2) 
where es is the baryon charge in units of e and kb is the baryon anomalous magnetic 
moment in natural units „, e r . 

An alternative representation of the electromagnetic form factors of the baryon 
B is the Sachs parametrization in terms of the electric charge G B and magnetic dipole 
Gm form factors. For the electric charge form factor the following relation holds with 
Fib and F 2B , 

Geb(Q 2 )=F 1b (Q 2 )-^F 2 b(Q 2 ). (2.3) 

As for the magnetic dipole form factor Gmb the natural definition is Gmb = 
Fib + F 2 b- At Q 2 — 0, Gmb(Q 2 ) defines the magnetic moment of the baryon B 
in natural units ( 2 m b )' l lB = Gmb(0) 2Mb ■ To compare magnetic moments of 
particles with different masses it is usual to express /is in terms [2n = 2Mn ; the 
nuclear magneton. In this case \xb = Gmb(0)jj^£in- Therefore, although the study 
of the baryon magnetic form factor can be done naturally using Gmb = Fib + F 2 b, 
as performed in a previous work [27], that \ib was defined in the natural units, it is 
more convenient to define Gmb in units of the nuclear magneton. Throughout this 
article we will use then 

Gmb(Q 2 ) - [F 1B (Q 2 ) + F 2B (Q 2 )] ^f- (2.4) 

M B 



2.2. In medium 

We consider now the octet baryon B in the nuclear medium with the effective 
mass M B . Assuming that the baryon is quasi- free in the nuclear medium, the 
electromagnetic current for the baryon B can be expressed as 

J'b = F* lB {Q 2 )r+F; B {Q 2 )^-^, (2.5) 
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where F* B and F£ B are respectively the Dirac and Pauli form factors in the nuclear 
medium. Again the in- medium spinors u* B (P + ) and u* B (P-) are suppressed. At 
Q 2 = 0, one has also 

F 1 * B (0)=e^, FS B (p) = K%, (2.6) 

where e* B is the electric charge in nuclear medium (the same as in the vacuum: 
cb = c* B ) and k* b is the anomalous magnetic moment in units of . 

As in the vacuum [see equations (|2.3|) - (|2.4|) ]. we define the electric charge and 
magnetic dipole form factors as 

G EB (Q 2 ) = F* lB {Q 2 ) - 1^2 F 2b(Q% (2-7) 

G* MB (Q 2 ) = [F* 1B {Q 2 ) + f; b {Q 2 )] (2.8) 

Note that the nucleon mass in vacuum (Mjv) is included in the definition of G* MB . As 
mentioned already we use this definition to make comparison easier with respect to 
the vacuum results. 

Because the effective nucleon mass is expected to be smaller than the mass in 
vacuum, Mfj < Mjv, G* MN (Q 2 ) is expected to increase and the magnetic moment is 
enhanced in magnitude (\n* N \ > |A*jv|)- 



3. Spectator quark model 

We describe now the octet baryon electromagnetic form factors in vacuum for a baryon 
B with mass Mb following [27] . In next section we describe the extension of the model 
to the nuclear medium. 

The electromagnetic interaction with a baryon B may be decomposed into 
the photon interaction with valence quarks, and with sea quarks (polarized quark- 
antiquark pairs or meson cloud). As the pion is the lightest meson the pion cloud 
is expected to give the most important contribution. Then, one can describe the 
electromagnetic interaction for a member B of the octet baryons using a current, 



J 0B t T J 7 £ 



(3.1) 



where J^ B stands for the electromagnetic interaction with the quark core without the 
pion cloud, and the remaining terms are the interaction with the intermediate pion- 
baryon (irB) states, as depicted in figure [TJ In particular, J£ represents the direct 
interaction with the pion [diagram (a)], and J^ B the interaction with the baryon while 
one pion is in the air [diagram (b)]. The factor Zb is a renormalization constant, which 
is common to each isomultiplet: nucleon (N), E, A, and 3. Zb is related with the 
derivative of the baryon self-energy 26J . 

We restrict the meson cloud dressing to the pion cloud, since the lightest meson 
is dominant as known from chiral perturbation theory. This is consistent with the 
studies of the octet baryon systems gillS|^|nil^|S3IMllS51|SBl|S7l|SS]. We note 
however, that kaon (K) cloud contributions may became more pronounced for systems 
with more strangeness, particularly when the pion cloud contributions are small. 

In the previous work [27| we presented a model for the valence quark and 
meson cloud contributions that were calibrated by lattice QCD and physical nucleon 
electromagnetic form factor data, as well as the octet magnetic moment data (A, 
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(a) (b) 

Figure 1. Electromagnetic interaction with the baryon B within the one-pion 
loop level through the intermediate baryon states B' . A diagram including a 
contact vertex ■yirBB' , as described in |26| . is not represented explicitly, since the 
isospin structure is the same as diagram (a). See 1261 for details. 



and S 0, ~). The model provided a good global description of the octet baryon 
electromagnetic form factor data (physical and lattice regimes) except for the neutron 
electric form factor in the low Q 2 region. As a consequence the model underestimated 
significantly the neutron electric charge square radius [—0.029 fm 2 to be compared 
with the experimental result —0.116 fm 2 ]. This was interpreted as an insufficient 
impact of the nucleon data in the low Q 2 region, particularly those for the neutron. 
Another limitation of the model was no explicit inclusion of the pion mass (m^) in 
the parametrization of the pion cloud dressing. Although the long-range falloff of the 
pion cloud effects in the position space with the distance r going like exp(— 2m 7r r), it 
can be simulated by multipole functions with appropriated cutoffs. In order to study 
the chiral behavior, it is vital to include the pion mass dependence explicitly. 

Therefore, we improve the model by adding two new features to the model of [27] : 

• We constrain the model parameters using also the experimental values of the 
proton and neutron electric and magnetic square radii, as well as E - electric 
square radius. 

• We redefine the pion cloud parametrization in order to reproduce the leading 
order chiral behavior for the form factors which depend on the pion mass. 

With these two additional constraints we can improve the form factors in the low 
Q 2 region consistently with the chiral behavior, in particular for the neutron electric 
form factor. We also provide a direct connection of the model with the chiral limit. 

Next, we describe how the valence quark and the pion cloud contributions are 
integrated in the model. The explicit parametrization for the pion cloud and the 
connection with chiral perturbation theory are presented later. 

3.1. Bare form factors 

In the covariant spectator quark model a baryon B is described as a system with 
an off-mass-shell quark, free to interact with photons, and two on-mass-shell quarks. 
Integrating over the two on-mass-shell quark momenta, we represent the quark pair as 
an on-mass-shell diquark with an effective mass uid, and the baryon as a quark-diquark 
system [20] [21] . This quark-diquark system is then described by a transition vertex 
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between the three-quark bound state and the quark-diquark state, that simulates the 
effect of confinement [21] ■ 

The simplest representation for a quark-diquark system with spin 1 /2 and positive 
parity is the S-wave configuration. As in [501 US 123 , we represent the wavefunction 
b for an octet baryon B with momentum P and the internal diquark momentum k. 
In the S-wave approximation the wavefunction is a combination of symmetric (\Ms)) 
and anti-symmetric (\Ma)) states in the quark pair (12), and an S-wave radial (scalar) 
wavefunction iPb(P, k). The explicit expressions can be found in |27j . 



3.1.1. Electromagnetic current. Taking into account that the wavefunction ^b is 
written in terms of the wavefunctions of a quark pair (12) and a single quark (3), 
one can write the electromagnetic current associated with the baryon B in impulse 
approximation (20j|3T] as: 

J ob = 3 E / * S ( P +' k)j^ B (P-,k), (3.2) 

P J k 

where j^ 1 is the quark current operator, P + (P_ ) is the final (initial) baryon momentum 
and k the momentum of the on-shell diquark. T = {s, A} labels the scalar diquark and 
the vectorial diquark polarization A = 0, ±. The factor 3 in equation Q3.2p takes into 
account the contributions for the current from the pairs (13) and (23), where each 
pair has an identical contribution with that of the pair (12). The polarization indices 
are suppressed for simplicity. The integral symbol stands for 

d 3 k 



2E D (2nf 



(3.3) 



where Eo = \J m 2 D + k 2 . 

Generally, the baryon electromagnetic current (|3.2j) can be expressed as 

in^ v n 

Kb = Fiob(Q 2 ) r + F 20B (Q 2 ) -^f , (3.4) 

where -Fios and i^os are respectively the valence quark contributions for the F\b{Q 2 ) 
and F2b(Q 2 ) form factors, defined by equation (|2.ip . To simplify our notations we 
introduce 

e 0B = F WB {Q 2 ), K 0B = F 20 b{Q 2 ), (3.5) 

where the tilde is included to remember that these are functions of Q 2 . To represent 
these quantities for Q 2 = 0, respectively the charge (Bob) and the anomalous magnetic 
moment (Rqb) we suppress the tildes as eoB and kqb- 



3.1.2. Quark electromagnetic current. The quark current operator j£ has a generic 
structure, 



< /! N icr'"'i h , 



tf=n^"^J +ia w' (3 - 6) 

where Afjv is the nucleon mass in vacuum and ji (i = 1,2) are SU(3) flavor operators 
acting on the third quark of the \Ma) or \Mg) state. In the first term 4q^ /q 2 is 
included for completeness, but does not contribute for elastic reactions. 

The quark current ji (i — 1,2) in equation (|3.6j) . can be decomposed as the sum 
of operators acting on the quark 3 in SU(3) flavor space [31] . 

ji = hfi+^O + + hfiO^s, (3-7) 
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Table 1. Mixed symmetric (jf) and antisymmetric (jf) coefficients for the octet 
baryons appearing in equations 1 13, 81 1 and 1 13. 9> . 



B 


it 




P 
n 


— fi— ) 
fi(/i+ + /»-) 


i(/ i+ + 3/<_) 
«(/»+ — 3/*-) 


A 




Ts(/»+ _ 4/ i0 ) 


£+ 
S° 
£- 


+ 3/<_ - 4/ i0 ) 

36( 2 /i+ - 8/jo) 

Ufi+ - 3/i- - 4/«,) 


£(/i+ + 3/<_) 
fi(/i+ — 3/,_) 


S° 


i(2/ i+ + 6/ i _-2/ i0 ) 
i(2/ i+ -6/ i _-2/ i0 ) 


-|/i0 
-f/iO 



where Ao = diag(l, 1, 0), A3 = diag(l,— 1,0) and A^ = diag(0, 0,— 2) are the flavor 
space operators. These operators act on the quark wavefunction in flavor space, 
q = ( u d s ) T . 

The functions fi±(Q 2 ) (i = 1,2) are normalized by fi n (0) = 1 (n = 0, ±), 
/ 2 ±(0) = K±, and /2o(0) = n s . The isoscalar (k+) and isovector anomalous 
magnetic moments are defined in terms of the u and d quark anomalous magnetic 
moments, k+ = 2k u — Kd and k_ = \k u + \n<i- 



3.1.3. Dirac and Pauli bare form factors. To take into account the effect of the states 
with mixed symmetry in the baryon B wavefunction (|Ma) and |Mg)) we sum over 
the quark flavors, using the coefficients 126} l27l [31] . 

■ ' ! - {M A \ji\M A ), (3.8) 



Jl 



jf = (M s \ji\M s ), 



(3.9) 



for i = 1,2. The expressions for jf and jf (i = 1, 2) are presented in table [T] -see [27] 
for more details. 

Using the coefficients defined by equations (|3.8|) and (|3.9[) we can write the 
spectator model form factors for the octet baryons characterized by a nucleon with 
mass Mjv and the baryon B with mass Mb as [2"T] 



with 



= Q 

117? 



cob 



hob 



and 



B(Q 2 ) 
B(Q 2 ) 



-i A 

2 Jl 



-J2 



13-t 



M 



-Jl 



21 + T 
11-3T 



2 1 + r 



1 + r M N 

Mb_ 2 
M N 



J 2 



1 



Ji 



B(Q 2 )= f tl> B (P+,k)MP-,k), 

Jk 



(3.10) 



(3.11) 



(3.12) 



the overlap integral between the initial and final scalar wave functions. 

The normalization of the wavefunction [20] 27| leads to B(0) = 1. Note that by 
construction the bare charge egs and the dressed charge es are the same, Cqb = £b- 
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We conclude then that the bare form factors S-qb and Hqb are determined by the 
fin(Q 2 ) (i — 1)2, n = 0,±) and i/j B (P,k). The details of those parametrization are 
shown next. 

We can also calculate the quark core contributions for the electric and magnetic 
form factors using the following expressions: 

Geob{Q 2 ) — e B - t R B , (3.13) 

G M0B (Q 2 ) = [e B +n B ]^f. (3.14) 
Mb 

In equation p. 141) the factor is included to be consistent with the definition of 
Gmb- 

Equations (|3~T0l> and (|3~TT|) can be used either in vacuum or in medium. In 
medium we may just replace the vacuum masses Mn and Mb respectively by the 
effective masses Mjy and M B . The same procedure should also be carried out in the 
coefficients jf and jf (i = 1,2), namely, the vector meson masses in vacuum should 
be replaced by those in medium. As for equations (|3. 13|) and Q3.14p the same rules are 
applied except that the factor in vacuum should be replaced by in medium 

(with the nucleon vacuum mass Mn) according to our convention for G* MB - 



3.1. 4- Parametrization of the quark current. To parameterize the quark current given 
by equation (|3. 61) . we adopt the structure inspired by the vector meson dominance 
(VMD) mechanism as in [2%ll3"T] . 



m 2 p MlQ 2 



fi± = \ + (1 - n2 + c ± n 



in 



Q 2 ^{Mi + Qiy 



f 1 o = A 9 + (l-A 9 )^- + c KQ2 



'{K + Q 2 ) 2 '' 

h± = K ± \d ± ^ + (1 Mh 



m 2 p + Q 2 v ±; A/2 + Q2 j ' 

where m p , and Mh are the masses respectively corresponding to the light vector 
meson (p meson), the cf> meson (associated with an ss state), and an effective heavy 
meson with mass Mh = 2Mn to represent the short-range phenomenology. We use the 
mass value m p for both isoscalar (corresponding u> meson) and isovector channels since 
m^j ~ m p . The coefficients cq,c± and dg,d± were determined in the previous studies 
for nucleon (model II) 20 and fi - [31]. The values are, respectively, c+ = 4.160, 
c- = 1.160, d+ = d- = -0.686, c = 4.427 and d = -1.860 [31]- The constant 
X q = 1.21 is obtained so as to reproduce correctly the quark number density in deep 
inelastic scattering [20]. 

The quark form factors parameterized by the VMD mechanism in equation (|3.15|) 
are particularly convenient to extend the model to other regimes besides the physical 
regime, because the quark current is written in terms of the vector meson and nucleon 
masses. In the previous work the model was extended successfully to the lattice 
QCD regime replacing these masses by those of the lattice regime [37]. Furthermore, 
the model was also applied to the lattice regime for the nucleon [3S], 7./V —> A 
reaction [3H [35] , and octet and decuplet baryons [2jj [31] . 
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3.1.5. Scalar wavef unctions. The scalar wavefunctions are given by [27] . 





N N 


(3.16) 


m D (Pi - 


^X»)(ft + X N )' 




N A 


(3.17) 


mnifix - 


HXa)(/3 3 + Xa)' 




iV E 


(3.18) 


mo (/3i - 


HxJ(ft + X s )' 






(3.19) 






m_D(/3i - 


^X S )(/?4+Xh)' 



where Nb {B = N, A, S, 5) are the normalization constants, and 

M B TOd 

Note that, except for the masses, the A and £ scalar wavefunctions are the same. The 
normalization constants Ng are determined by 

|^(P,fc)| 2 = l, (3.21) 

k 

where P = (Mb, 0, 0, 0) is the baryon four-momentum at its rest frame. 

In equations p,16|) - (|3.19|) the parameters /?,; (i = 1,..,4) define the momentum 
range in units of mp. The parameter f3\ is associated with the long-range scale 
(low- momentum range) that is common to all the octet baryon members. As for the 
remaining parameters /32, Pz, Pi > they are associated with the shorter range scale 
(larger momentum range). Namely, /?2 defines the short-range scale for the systems 
with only light quarks u and d, /?3 defines the short-range scale for the systems with 
one strange quark, and (3^ defines the scale for the systems with two strange quarks. 
As the strange quarks are heavier than the u and d quarks, and therefore more confined 
in the space, we expect that /?2 > fiz > Pa- The parameters /?,; (i = 1, .., 4) as well as 
k u , Kd will be fixed later. 



3.2. Pion cloud dressing 

We discuss here the pion cloud contributions for the electromagnetic current and form 
factors represented by the diagrams in figure [U Following |26| , we assume the pion as 
the dominant meson excitation to be included in the octet baryon form factors. Then, 
the meson cloud contributions for the octet baryon electromagnetic form factors can 
be described in terms of 6 independent functions of Q 2 , related to the pion-baryon 
Feynman integral as will be described next. 



3.2.1. Pion cloud electromagnetic currents. The pion cloud corrections, namely the 
coupling of the photon to the pion J£, and the coupling to the intermediate baryons 
J^ B , defined by equation (|3.1[) . can be written [26l [27] 

J! = (BiY + B,'^^-) G,b, (3.22) 
V + (3.23) 
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In the above, Bi, Ci and Di (i — 1,2) are functions of Q 2 and G^b, G s b, G k b are 
coefficients that depend on the baryon species (B = N, E, A, 2). We assume that the 
functions Bi, Cj and Di are only weakly dependent on the baryon masses, and the 
same for all the octet baryons as in [26 . That allows a description of the pion cloud 
dressing with a reduced number of coefficients. We write Bi, Ci, and Di to represent 
respectively the functions Bi, d and Di at Q 2 = 0. 

The coefficients G^b, G c b, G k b include the dependence on the pion-baryon 
coupling constants. According to SU(3) symmetry [59l [60] the coupling constant of 
the pion (%) and baryons {B and B'), g-xBB', can be represented in terms of the ratio, 
a = p ® D and a global coupling constant g = g^NN, the nNN coupling constant. We 
can then express G n B, G c b, G k b in terms of the parameter a and a global factor 
g 2 . For convenience we absorb the factor g 2 in the functions Bi, Ci and Di (i = 1, 2) 
and represent the effect of the coupling in terms of 4 independent constants [26l [27] 
associated with the octet baryon specie^], 

Pn = 1, (3.24) 

Pa = \a 2 , (3.25) 

/3 E =4(l-a) 2 , (3.26) 
f3 3 = (1 - 2a) 2 . (3.27) 

These constants encapsulate the effect of the coupling constants. The explicit 
dependence of G^b, G e s and G k b on the constants given by equations (|3.24[) - (|3.27j) 
and on the bare form factors e^s and kqb was derived in J26[[27]. In the following we 
use the results of [27] . 

From the equations above, we get (3n = 1, /?a = 0.48, /3s = 0.64 and /3~ = 0.04 
with a = 0.6, determined in combination with an SU(6) quark model. It is therefore 
expected that the pion cloud contributions are small for the S system. In this case 
the kaon cloud contribution may be more significant. 

3.2.2. Dressed form factors. The octet baryon dressed form factors associated with 
the current (|3.ip . are obtained by the contributions from the quark core given by 
equations f|3 . 10[) - (|3 . 1 1[) . and the pion cloud dressing via equations (|3.22p - (|3.23j) : 



J£ = Z B 



cob + GttbBi + G e sCi + G k bD\ 

KOB + GhbB2 + G e bC2 + G K sZ?2 



7 M + 



2Mb 

(3.28) 

Using the above expressions and the definition of the form factors (12.1[) . we can write 
down the final results for the form factors Fib and Fib ■ 

Fx p = Z N |e 0p + 2f3 N Bi + /3n(c 0p + 2e 0n )Ci + /3 N (k 0p + 2k 0n )Di j , (3. 
Fin = Z N |e „ — 2/3 N Bi + f3 N (2e 0p + e 0n )Ci + /3 N (2k 0p + k 0n )L>i j , (3. 

Fia = Z A |e A + /3A(e s+ + e 0E o + e os - )C\ + /3a(«os+ + ^os° + «os- )D\ 

| The coefficient /3jv was not considered explicitly in the previous works 26 27 , where f$pj = 1, but 
it is included here for completeness and to clarify the extension of the model to the in-medium case. 
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(3.31) 



= {g os+ + (/3 E + /3 A )Bi + [/Ss(Sos+ + S 0S o) + ^aSoa] Ci 

+ [/3s(^os+ +koe«) +/3a^oa]£»i} , (3.32) 

F 1E o = Z s |e os o + [^ s (e E+ + e s-) + /SaSoa] Ci 

+ [/3e(^os+ +kos-)+/?akoa] j Di}, (3.33) 

^ie- - {Sos- - (/3e + /8a)Si + [/3e(S 0S " + eos-) + /?aSoa] Ci 

+ [/3e(^oso + ^os- ) + /3a^oa] £»i } , (3.34) 

F ls o = {g OH o + 2/3 H Bi + ^ S (e 0H o + 2e 0H -)Ci 

+/fe(£ 0H o +2£ 0H -) J Di}, (3.35) 
F 1H - - Z s {g 0H - - + /fe(2g OH o + e 0H -)Ci 

+^ h (2k 0H o + k oh - )D 1 } , (3.36) 



F 2p — Z 



n |«o P + 2/3 N B 2 + /3 N (eo P + 2e 0n )C 2 + /3jv(«o p + 2k 0n )D 2 } , (3.37) 
^2n = {«0n - 2/3jv£j 2 + M2e 0p + e 0n )C 2 + f3 N (2k 0p + K 0n )D 2 } , (3.38) 

F 2 \ = Z A |k a + A\(e £+ + e E" + e £-)C 2 + /3a(k os + + k e° + Koe-)-^} , 

(3.39) 

^2S+ = ^e {«os+ + (/?£ + ^a)-B2 + [/fe(e os + + e os o) + /3 A e 0A \ C 2 

+ [/3e («oe+ + «os° ) + j9a«oa] D 2 } , (3.40) 

F 2S o = Zy, {k £° + [^s(e s+ + e £-) + /3aCoa] C 2 

+ [/3e (k e+ + kos- ) + /3a«oa] D 2 } , (3.41) 

F 2E - = Z s |koe- - C#£ + /3\)B 2 + [/fe(e 0E o + e s- ) + Ay^a] C 2 

+ [/3s (ft E" + «oe- ) + PaHoa] D 2 } , (3.42) 

F 2S o = ^e {koh« + 2/3 s S 2 + /3h(Soho + 2e 0E -)C 2 

+ f3 s (k 0S «+2k 0E -)D 2 } , (3.43) 

F 2H - = Z E {k oh - - 2^ H S 2 + /3 H (2e OH o + e 0H -)C 2 

+^ h (2k oh o + H os - )D 2 } . (3.44) 
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The normalization constants are 
r n -1 



Zn — 


1 - 






Z A = 


H 


- WkBx 


-l 


Zv = 


1 - 






z s = 


H 


- 3/3 S -Bi 


-1 



(3.45) 



and we refer to [251 [S7J for more details. 

Using the above expressions, one can calculate the electric and magnetic form 
factors both in vacuum, based on equations (|2.3I) and (|2.4I) . and in medium based on 
equations (|2~71) - (|2~5)) . 

From the discussions in the previous sections, the baryon form factors Fi B 
(i = 1, 2) may be decomposed into 

Fw(Q 2 ) = Z B [F MB (Q 2 ) + 5F iB (Q 2 )] , (3.46) 

where -Fios and i^os are defined by equation (13.41) . and the corresponding pion cloud 
contributions SF 1B and 5F 2B are given in equations (I3.29[) - (I3.44[) . It is natural to 
regard Z B F^, B as representing the valence quark effects and Z B 5Fi B those of the pion 
cloud. 

The same decomposition can be applied for the electric and magnetic form factors: 
G EB (Q 2 ) = Z B [G E0B (Q 2 ) + 5Geb{Q 2 )} , 

G MB (Q 2 ) = Z B [G M ob(Q 2 ) + SG MB (Q 2 )] , (3.47) 

where Geob,Gmob are defined by equations (|3.13l) and (I3.14[) . and SGeb — SFi B — 
tSF 2B and 5Gmb — [SFib + SF 2B ] ^f-. In this case Z b SGeb, and Z b SGmb reflect 
the dressing of the pion cloud. To estimate the pion cloud contributions, we compare 
the full result, Geb or Gmb, with the total contributions of the valence quark core, 
ZbGeob or Z b Gmob- The difference is the pion cloud contributions, Z b SGeb or 
Z b 5Gmb- 



4. Including chiral symmetry in the pion cloud parametrization 

The effect of chiral symmetry in the electromagnetic structure of the octet baryons 
can be analyzed by studying the dependence of the form factors on m w and the 
corresponding behavior in the chiral limit (m, —>()). In the chiral limit the pion 
cloud extends to infinity, leading to the divergence of the nucleon radius. 

In this section we start reviewing the main features of the chiral behavior for 
the nucleon radii. At this stage we do not attempt to describe the nucleon magnetic 
moments including the pion mass dependence, since the dependence is milder than 
that for the nucleon radii [BT1 IEH Next, we derive the general expressions 

for the nucleon radii in the present model. Finally, we present the newly updated 
parametrization for the pion cloud contributions, and describe how the chiral behavior 
is implemented in the model. 
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Table 2. Nucleon radii. We adopt the average values from the PDG 1671 more 
recent results. All values in fm 2 . The Pauli isovector square radii normalized is 

^2^- = 0.418 ± 0.012 fm 2 . 

"71 („v \i i„v \i 
'Jig ' Mp ' En ' Mn V 1 ) V 2 / 

0.7634±0.0140 0.6983±0.0109 -0.1161±0.0022 0.7430±0.0134 0.0634±0.0142 1.5509±0.0430 



4-1. Nucleon radii in the chiral limit 

The leading order effects of chiral symmetry can be better observed in the nucleon 
isovector form factor [64] defined by, 

F^(Q 2 ) = F lp (Q 2 )-F ln (Q 2 ), (4.1) 

F% (Q 2 ) = F 2p [Q 2 ) - F 2n [Q 2 ). (4.2) 

Using these decompositions we can define the isovector Dirac and Pauli square radii: 

dFY 



{rXf = -6 



dQ 2 



(4.3) 

Q 2 =0 



( r 2 ) = -6 



dQ 2 



(4.4) 

Q 2 =0 

Note that in the definitions we do not normalize the respective radius at Q 2 = 
as usually done. To compare (rY) 2 with the experimental value the result must be 
divided by Ky = k p — K n ~ 3.7 (isovector anomalous magnetic moment). 

We now discuss the expected result in the small pion mass limit. According to 
XPT |65j|66] the isovector square radii (rY) 2 and (r^) 2 can be expressed as 

(rf) 2 = -— logm^+^i, (4.5) 

(rl) 2 = + ^ + A 2 , (4.6) 

where and M are the pion and nucleon masses, respectively, and 

a = 8ir 2 F 2 , (4.7) 

ai = 5g 2 A + 1, (4.8) 

a 2 = irg 2 A , (4.9) 

with F,r = 93 MeV and g A — 1.27. In the above, Ai and A 2 represent constant terms 
and higher powers of m 2 , that can be expressed by constants at the physical point. 

We take a pragmatic approach and fix these constant values using the 
experimental values of (rY) 2 and (rY) 2 - The results obtained from the average valued 
of the PDG results [57] are presented in tabled From the table we find A\ ~ 12.05 
fm 2 and A 2 ~ —38.20 fm 2 . These results are obtained assuming the physical value of 
m-n and the experimental results. 

§ For the neutron radii r% and r'^ ln we use the results suggested by PDG 67 . For the proton radii 
there is still controversy. For rj v[p we take the arithmetic average of the 3 results listed. As for r|, 
we simply take the arithmetic average of the first 12 results from the list (since 2000), irrespective of 
the different types of the determinations, electronic, muonic, or others. 
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4-2. Isovector form factors and respective radii 

We discuss now the isovector form factors given by the present model. 
Using equations ([3T2l?|) - ([3~5(Ij) and {537]) - (g35) we can write 

FY{Q 2 ) = Z N {(S 0p - e 0n ) + 4Bi - (i Qp - g 0n )(7i - (k 0p - k „)^4.|q) 

F 2(Q 2 ) = z n {(^o P - H 0n ) + 4B 2 - (e 0p - e 0n )C 2 ~ {k 0p ~ £ n)(iP-i|) 

From these equations we conclude that the nucleon isovector form factors are given 
by the difference between the proton and neutron bare form factors (Son or kojv) and 
the pion cloud contributions. 

The results for the isovector square radii are given by 



where 



(rf) 2 = 

i?i = 6 
i? 2 = 6 



J N 



z 



N 



-24 



-24 



dB 1 

dB 2 
dQ 2 



Ri 



dCx 

w 

dC 2 



+ (1 - Ci){r 2 Wp - r\ Qn ) + 6(n 0p - n Qn ) 



dD x 



dQ 2 



(4.12) 
(4.13) 

(4.14) 



(r 



2 

10p 



dQ 2 

■ 6(k 0p - KOn) 



-r 2 Wn )C 2 + (l-D 2 )(r\ 
dD 2 



2 

20p 



dQ 2 



In the above expressions | n stands for the derivative at Q 2 — 0, and 



' iOB 



-G 



dF, 



iOB 



dQ 2 



(<=1,2), 



(4.15) 



(4.16) 



3 2 =0 



represent the bare radii. 

The expressions ()4.12p and ()4.13p are still general. We now discuss the constraints 
of chiral perturbation theory for the model. 

4-3. Pion cloud "par ametrization 

We consider the following parametrizations for the functions B\ and B 2 : 



B 1 = B 1 



£>2 — B 2 



A? 



Aj + Q 2 

A 2 



1 



1 



Z N B X V24a 



1 Cki 



log m w + b[ I Q 



6 r 



1 



ZnB 2 



1 a 2 M 
24 ofo m v 



b' 



(4.17) 
(4.18) 



Here B\ and B 2 are constants given respectively by Bi(0) and B 2 (0), and Ai, A 2 arc 
two cutoffs to be fixed by a fit to the data. 

As for b[ and b' 2 they are two additional parameters that will be fixed by the 
experimental results for the nucleon isovector square radii, equations (I4.5[) and (I4.6[) . 
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Inserting the above expressions into equations (|4.12l) and (|4.13p we can write 
{r vf = _2l i ogm7r + Zn (-24A- - 120§ + ife ) , 



a m n { Z N 



(4.19) 

v^a _ ^ -M ^ 7 _ /_ ,A _ 144 |l _|_ . 

(4.20) 

Note that the second term in each equation above should be identified with A\ and 
A2 in equations (|4.5I) and (|4.6I) . respectively. 
As for the remaining functions we use 

Ci=Bi[jr^s) , (4-21) 
^ 1 = ^(A^W' (422) 



a: 



Again C2 and D2 are the constants given by the value at Q 2 = for the respective 
functions. D[ is a new constant defined by D[ = Note that Ci(0) = B x , 

a constraint required by the conditions -Fip(O) = 1 and Fi„(0) = (nucleon 
charges) [25] [27J. The definition of D\ that vanishes at Q 2 = is also motivated 
by the nucleon charge conditions. The parametrizations used for C\, Di, C 2 and D 2 
are the same as those presented in [27]. The leading order chiral effects in the form 
factors in the present parametrization come exclusively from B\ and B 2 . We choose 
to use the same cutoffs in Ci,Di as those used in Bi (Ai and A2, for the Dirac and 
Pauli form factors) in order to reduce the number of parameters in the model [27]. 
Finally, we can now write down, 



/?j = - 12^ + (l-Ci)(r?op-r? 0n )+6(/s Qp -/so„)^ ) (4.25 i 



R2= -18§ + (l-^2)(r 2 2 o p -r 2 2 J 



Do 



- ( r io P - r wn)C2 - 18(k 0p - «on)-r|- (4-26) 

A 2 

The values of the coefficients b[ and b' 2 can now be determined comparing equations 
P~T9l and P~20)) with equations (g^J) and (j4~cll) . 

The choice of the powers included in the pion cloud functions is phenomenological 
and motivated by the expected falloff of the quark-antiquark contributions in the large 
Q 2 limit |68] as well as the magnitude of the pion cloud contributions estimated for 
the jN — > A reaction [24] [33] [34] . With the present parametrization the pion cloud 
contributions fall off by a factor 1/Q 4 faster than the falloff of the valence quark 
contributions. 
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5. Extension of the spectator quark model to the in-medium regime 

We now discuss the extension of the model for the in-medium regime. In general 
we consider the modifications of the model due to the in-medium environment. As 
for the valence quark core part, the in-medium hadron masses appearing in the 
model, Mb (B = iV,A,S,H), m p , and Mh, will be respectively denoted by M B 
(B = N,A,E,E), to*, toJ and M£ (given by 2M^). On the other hand, for the 
pion cloud effects we consider the modifications of the pion-baryon couplings in the 
in- medium regime as will be explained next. In this study we do not include any final 
state interactions nor meson exchange current, as mentioned in introduction. 

5.1. In-medium regime: quark core 

Although lattice QCD simulation has been very rapidly developing recently, it is still 
very difficult to study the properties of hadrons in finite nucleon (baryon) densities near 
and higher than the normal nuclear matter densities. Thus, we need to resort to some 
phenomenological models which have proven successful in studying nuclear phenomena 
and nuclear processes based on the quark degrees of freedom. For this purpose, we use 
the quark- meson coupling (QMC) model [39j[40], which has been successfully applied 
to study the properties of baryons and mesons in nuclei and nuclear medium. 

We note that, combined with the cloudy bag model (CBM) [55], QMC had 
indeed predicted [TU] an in-medium modification of the bound proton electromagnetic 
form factors which turned out to be consistent with the experimentally observed 
modification [U [5j [6j [7]. Below, we use a different model, the covariant spectator 
quark model [20] [23l [24] , which was successfully applied to study the octet baryon 
electromagnetic form factors [27] utilizing the lattice QCD simulation data. By the 
use of different models, we hope to shed light on the mechanism of the in-medium 
modification of the bound proton electromagnetic form factors. Furthermore, we 
predict in-medium electromagnetic form factors of all members of the octet baryons. 

The current j^, given by equation (|3.6p . is also characterized by the corresponding 
in- medium masses in terms of the two components j\ and j'2, which are represented 
based on the VMD parametrization in equations (|3.15|) . 

In the quark current (|3.6[) , we replace the coefficient of the Pauli form factor 
1/(2Mjv) by 1/(2M^) in the in-medium regime. As for the quark form factors, we use 
equation (|3.15|) with the meson masses replaced by the respective in-medium masses 
similarly to the lattice regime studies [27l [3TJ |3H [35] . That is, we replace m p and 
by the in- medium p mass to*, and </> mass to^, and the effective heavy meson mass of 
M h = 2M N by 2M^. 

As for the wavefunctions "J b , Mb is replaced by M B in medium. This applies for 
the radial (scalar) wavefunctions (|3.16l) - (|3.19[) . As explained the scalar wavefunctions 
for the octet baryons are represented in terms of four independent momentum range 
parameters (i = 1,..,4). We assume that these parameters are independent of 
the baryon masses in vacuum, and therefore also independent of the in-medium 
baryon masses. There is no need to modify the diquark mass in in medium, since 
the electromagnetic form factors are independent of it [20j [27] . For the in-medium 
regime the pion mass to*, which is relevant for the present study, it was estimated 
in [72] that to* ~ to^ at normal nuclear matter density (0.15 fm -3 ). Thus, we use 
the vacuum value, = 138.0 MeV for the densities considered in this study. The 
in- medium hadron masses mostly calculated in QMC, and the relevant values for the 



Octet baryon electromagnetic form factors in nuclear medium 



17 



Table 3. Hadron masses (in MeV) necessary for the in-medium regime of the 
model with po = 0.15 fm — 3 . For the vacuum case (p = 0) we take m p as the 
average of p and u masses as originally used in the model 1201 127| . 





P = 


p = 0.5,00 


P = Po 


M N 


939.0 


831.3 


754.5 


M A 


1116.0 


1043.9 


992.7 


M E 


1192.0 


1121.4 


1070.4 


M s 


1318.0 


1282.2 


1256.7 


m p 


779.0 


706.1 


653.7 




1019.5 


1019.1 


1018.9 




138.0 


138.0 


138.0 



calculation, are listed in table [3J 

Using the model extended to the in-medium regime, namely, using the quark 
currents and baryon wavefunctions for the in-medium regime, we can calculate the 
form factors G* EB and G* MB in a nuclear medium by the expressions given in section 
13.11 for a model with no pion cloud dressing with the in-medium masses calculated in 
QMC for a given nuclear density. 

5.2. In-medium regime: pion cloud 

To extend the pion cloud effects of the model to the in-medium regime, we need other 
ingredients. Although we do not have a good control, we estimate the effects of the 
pion cloud in the nuclear medium. Thus, the estimates of the pion cloud effects in the 
nuclear medium presented below, should be taken in caution. Main task for this is to 
calculate the modifications of the baryon-pion coupling constants in medium, g* BB , 
(B, B 1 = N, A, S, H). [As before, we denote the properties in the nuclear medium by 
asterisk *.] For the moment, we consider below the ir BB diagonal case, and omit the 
7rAE coupling which appear in the pion cloud effects in the A and S cases containing 
the g^AS coupling from the discussion, but the procedure can be extended to also 
for the 7rAS case. We assume that the pion cloud parametrization functions given 
by equations (|4.17[) - (I4.18I) and (|4.21[) - (l4.24j) defined in vacuum, are unmodified in the 
medium, since the pion mass in vacuum is used. 

Our estimates of the in- medium couplings, g^ BB , relative to those in vacuum 
g-nB b [69| I70j , rely on the Goldberger-Treimann relation [71j . The in-medium to the 
free coupling constant ratio may be expressed by, 

9i BB = ([A (9a*\ rm\ 

9,bb \f*J \g% ) \M B )' 

where /„■, g^ and M B are respectively the pion decay constant, axial coupling constant 
of the baryon B and its mass, and the corresponding quantities in nuclear matter with 
*. First for /*, it was estimated in [75] and we use the f^ in nuclear matter as /* 
above. On the other hand, g^*/9A was estimated in QMC [73], and M B /M B as in 
tabled For g*^, the relevant quantities different from the diagonal cases are M^/Ma 
or M^/Afs- To estimate the maximally modified case, we use M^/Ma, although the 
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p = p = 0.5p P = Po 



91nn/9hNN 



1 0.921 0.899 

1 0.973 0.996 

1 0.977 1.004 

1 1.012 1.067 



Table 4. Modification of the txBB' 
po = 0.15 mi" 3 . 



coupling constants in nuclear matter with 



difference is less than 1%. The modification of the coupling constants for the densities 
p = 0.5 po and po are summarized in table 01 In table 0] the in-medium coupling 
constants g^BB' YVA equation (|5.ip either decrease, or remain close to the vacuum 
values. In QMC it is expected that the values slightly decrease in medium, since the 
coupling is the same as that of the weak axial coupling constant gA, namely 7^75 to 
the Dirac spinor, and g^* in nuclear matter decreases [73| . Some unexpected behavior 
for g* SE and <?* ss are due to the large decrease of /* estimated in [72] . The decreasing 
rate of /* in medium overcomes slightly that of the g^*M^ and g^*M^, and thus 
3tt££ an d ff^HH increase slightly. Thus, the trend of changing in values for g* B B' 
in nuclear matter, is consistent with that expected from QMC, although the latter 
generally does not contain the pion cloud effects, and it is this pion (/*) which leads 
to slightly unexpected density dependence for g* ss and g* ss - 

Next, we comment on the effects of the baryon mass modifications in the 
intermediate baryon propagator in the nuclear medium. In the pion cloud dressing of 
the current shown in figure [I] in vacuum, the intermediate states baryons belong to the 
same isomultiplet except for the A and £ cases. Therefore the intermediate baryon 
masses are the same or at most differ 1% (A and S). In-medium mass modifications 
of the baryons B' apply in the same way as those for the initial and final baryons, 
and thus the mass difference between the baryons B and B' in the nuclear medium 
is the same as those in vacuum, namely, [Mb 1 — Mb] — [M B , — M B ], and we can 
approximately have the same pion cloud effects due to the modifications of the baryon 
masses. Thus, the modification of the pion cloud effects in the nuclear medium arises 
entirely from the modification of the pion-baryon coupling constants, g^BB' > m ^ ne 
present approach. 

Finally, we estimate the in-medium modifications for the pion cloud dressing by 
replacing the constants those in vacuum /3b (B = N, A, S, S) given by equations (|3.25p ~ 
(I3.27|) . by the new constants j3* B to be presented next. We keep the value a — 0.6, the 
same as that in the vacuum, and use the in- medium coupling constants g^BB'i an< ^ 



get f3* B for B = N,A,E,E, 




(5.2) 



(5.3) 



(5.4) 



(5.5) 
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Note that according to the relation between the constants fig and the normalization 
constants Zb [see equations (|3.45l) ] the renormalization constants will also be modified 
to Z* B (replacing j3 B -> j3* B ). 

6. Results 

According to the new constraints imposed on the model (chiral symmetry and fit to 
the experimental square radii) we readjust the parameters of the model and carry out 
a new calibration that differs from the one presented in [27] . 

We divide the presentation of our results in two parts. First, we will present the 
results of the new calibration for the vacuum, mainly concentrating on the nucleon 
form factors. Next, we will present the results for the baryon octet electromagnetic 
form factors in the nuclear medium, showing also each contribution from the valence 
quarks and pion cloud. 

6.1. Octet baryon electromagnetic form factors in vacuum 

We present here the results of the global fit in vacuum, and compare with the 
previous model in [27j . The values of the valence quark parameters including the 
momentum-range scales of the octet baryon wavefunctions and the u {k u ) and d {nd) 
quark anomalous magnetic moments are presented in table [S] As for the parameters 
associated with the pion cloud, they are given in table [51 

We fit the parameters of the model to the nucleon data [74 ] |75] ITS ] l77 l [78j 179 ] [80j 
EH HU IH31 IHH HZ], the octet lattice QCD data, the octet magnetic moments 

(A, S + '~, S°'~) [57], the nucleon electric and magnetic radii and also the E~ electric 
radius. Following [57J we use some constraints in the fit to the octet lattice and 
physical data. The lattice data considered are the data with m w = 351, 591 MeV from 
[88] for n,p, E 1 * 1 , S 0,_ . For these pion mass values, the pion cloud effects arc expected 
to be small. 

To perform a fit we have taken care of the following points in order to achieve a 
fair description of both the physical and lattice data: 

(i) The impact from the nucleon physical form factor data is doubled compared with 
the octet lattice form factor data. 

(ii) Statistical errors are doubled for the neutral particles n and S° to take into 
account the possible systematic errors in the lattice. 

(iii) We double the error bars for the physical magnetic moment data to avoid 
dominating \ 2 due to the extremely accurate measurements in comparison with 
the nucleon physical form factor data. 

(iv) The impact from the nucleon radii data are doubled to be consistent with the 
nucleon form factor data. 

(v) The experimental error for the S~ electric square radius is reduced in our fit 
(from 0.15 fm 2 to 0.015 fm 2 ) in order to avoid the dominance of the nucleon 
square radii data in the fit and also to get the same order of contribution from 
the pion cloud as that for the proton. 

Below we explain in detail the fitting procedure described above. 

The first condition, the reinforcement of the impact of physical data in the fit 
is included since the number of physical data points (202) is inferior to the number 



Octet baryon electromagnetic form factors in nuclear medium 



20 



of lattice data points (272) as discussed in [37]. The reduction of the impact of the 
neutral particle lattice data was also discussed in detail in [37] . As for the magnetic 
moment data (magnetic form factors for Q 2 = 0) the fitting condition is included 
to avoid the excessive dominance of the Q 2 = data. The first 3 conditions follow 
the procedure used in the previous work where the nucleon physical data and octet 
magnetic moments and the octet lattice data were used to calibrate the model |27j . 
We discuss now the inclusion of the data related with the electric and magnetic radii 
that explain the last two conditions. 

As discussed, the available experimental information about the octet radii is 
restricted to the nucleon (p and n) electric and magnetic radii and the E _ electric 
radius. The nucleon square radii are given in table [2] For the E~ electric square 
radius the experimental value is r^ s _ = 0.61 ± 0.15 fm 2 89 . In the \ 2 calculation 
we use a factor 2 for the radii data to be consistent with the physical and lattice 
data in the \ 2 evaluation (impact 2, the same as the nucleon form factor data). The 
inclusion of the octet radii data is important to calibrate the low Q 2 regime of the 
model. In addition, the radii data constrain the degrees of freedom related to the pion 
cloud dressing, and reflect the effect of the chiral symmetry. Note however the large 
error associated with the E~ electric square radius. Although the nucleon radii data 
restrict the possible values to a small interval, the E~ electric square radius has a 
large interval of variation. 

The minimization of the x 2 using the result for the E _ electric square radius 
(with the error 0.15 fm 2 ) leads to a very large E~ electric square radius compared with 
that of the proton, as a consequence of much larger contribution from the pion cloud 
than that for the proton. A large pion cloud contribution for the E~ electric square 
radius contradicts what it is expected from SU (3) symmetry, where the contribution 
is expected to be of comparable with that for the proton. Also, estimates of the pion 
cloud contribution from |48] for the electric square radii predict similar contributions 
for the proton and E _ . In order to impose a small pion cloud contribution to the E _ 
electric square radii with the amount close the the proton case, we reduce the error 
of E~ electric square radius in the \ 2 calculation to 0.015 fm 2 (instead of 0.15 fm 2 ). 
We consider then, for the E~ electric square radius an error bar comparable with the 
error bars from the nucleon radii. With the choice 0.015 fm 2 , the contribution of the 
E~ electric radius to x 2 is of the same order as the other contributions (particularly 
the proton magnetic square radius), and the calibration of the model is not dominated 
by a particular observable. An additional motivation to use a standard deviation for 
the E~ electric square radius with a magnitude ss 0.02 fm 2 is that we can achieve 
a simultaneous good description of the nucleon form factor data (x 2 per data point 
2) and lattice data (\ 2 per data point ss 6), as in the previous work, where the 
experimental information about the baryon radii was not taken into consideration |27) . 
We note that this is an ad hoc procedure, but in the absence of additional information 
about the quark core effect (or pion cloud) in the physical regime, it is the simplest 
way to make a realistic calibration of the present model. Later we will discuss the 
sensitivity of the final fit to the values of the T,^ 1 electric radii. 

From the parameters associated with the valence quark contributions in table 
13. 1| we note that 0i < 02,03,04 as expected from the interpretation of fix, as the 
3-quark long range parameter according to equations (I3.16[) - (I3.19[) . Also the order of 
the momentum range scales 02 > 03 > 04 are consistent with the fact that the system 
with only light quarks (parameter 02) is more spread than a system with one strange 
quark (parameter 3 ), and that it is less compact than a system with two strange 
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Pi 


Pi 


Pa 


Pi 


0.0532 


0.809 


0.603 


0.381 








K s 




1.711 


1.987 


1.462 



Table 5. Parameters associated with the valence quarks determined by the 
fit. The value of k s was determined in |31| in the study of the baryon decuplet. 
The analytical form of the scalar wavefunctions, depending on fti (i = 1, .., 4) are 
presented in section [XT] 



Bi,B 2 


c 2 


D[,D 2 


b[,b> 2 


Ai,A 2 (GeV) 


0.0510 




-0.148 


1.036 


0.786 


0.216 


0.00286 


0.0821 


-1.987 


1.132 



Table 6. Parameters associated with the pion cloud dressing determined by the 
fit. See the paramctrization of the functions Bi,Ci and Di [i = 1,2) in section 
13.21 The values of b^b^ are calculated from the equations Il4.5l l-l l4.6l l. 

quarks (parameter ^4) in the position space. The same trend was observed in |27j . 

As for the values of the pion cloud parameters the more significant difference is 
the increase of Bi, which implies an increase of the pion cloud contribution for Ge, 
for the all members of the octet baryons, near Q 2 — 0. Taking the nucleon case as 
an example, Zm = 1/(1 + 3£?i) ~ 0.87 in the vacuum, means that the pion cloud 
contribution for the proton charge is about 13%. This is still rather small compared 
with a model such as CBM, yet in comparison with the previous work, the low Q 2 
behavior of the functions B\ and B 2 is modified. (See section |4~3"|) . This modification 
changes significantly the behavior of the pion cloud parametrization. Nevertheless, 
the new parameter values are close to the previous ones except for C 2 |27j . which is 
now much smaller (reduction in the effect of the Dirac form factors due to the process 
shown in figure QJb) to the magnetic form factors). 

Comparing the quality of the present fit with the previous one, we have obtained 
a slightly worse description of the lattice QCD data [x 2 per data point of 6.0 to be 
compared with 5.0], and also a less accurate description of the nucleon data \\ 2 per 
data point of 1.99 to be compared with 1.93]. In detail we have now per data point, 

X 2 (Ge p ) = 1.89, X*(G M p) = 1-69, 

X 2 (G En ) = 1.78, X 2 {G Mn ) = 2.41, (6.1) 

to be compared with the previous values of X 2 {Ge p ) — 1-60, X 2 (Gm p ) = 1-87, 
X 2 {Geu) — 1-86 and x 2 (Gmu) = 2.27 (27]. In simple words we have improved the 
description of Geu but we have lost some precision in the description of the other 
form factors, particularly for Gmu- The neutron magnetic form factor is difficult to 
describe as a consequence of the high accuracy of recent data, that disagree with the 
previous sets (see details in [27]). The loss of precision in some nucleon physical data 
is compensated by the quality of the description of the nucleon and £~ square radii, 
where the constraints from chiral symmetry are imposed. 

The main electromagnetic properties of the octet baryons in vacuum, such as 
magnetic moments (in /ijv units), the electric and magnetic square radii are presented 
in Table [7] and compared with the experimental results. 
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-0.689 


-0.65(3) 
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Table 7. Electromagnetic proprieties of octet baryons in vacuum. See Table [2] 
for the description of the nuclcon radii data. Additional data are from [89) (S _ 
electric square radius) and 67 (magnetic moments). 



The square radii are denned according to 
2 6 oIGeb 



' EB 



G EB (0) dQ 2 
6 dG M B 



(6.2) 

Q 2 =0 



' MB 



Gmb(0) dQ 2 



(6.3) 

Note that the square radii are normalized by the value of the form factor (Geb or 
Gmb) at Q 2 = 0, the usual definition. For neutral particles [Geb(0) =0], we use the 
same definition with Geb(0) — > 1. 

In table [8] we present also the decomposition of the octet baryon square radii 
into valence or bare core (b) and pion cloud (tt) contributions. These are defined for 
X = E,M, by 



2 \ _ ry 6 dGxoB 



( r XB)b =~ Z B 



G XB (0) dQ 2 

, 2 >> 7_ 6 d(SGxB) 
vxb) 



(6.4) 

Q 2 =0 



(6.5) 

2 =o 



"Gxb(0) dQ 2 

based on the decomposition (|3.47[) . Note that the components {r 2 XB )b and (r^ B ) OT are 
normalized by the total form factor Gxb(0). 

In the fitting process we can conclude that the model is very sensitive to the 
neutron data. This can be a consequence of the impact of the pion cloud in the 
neutron form factors. [This can also be true for other neutral particles, but there is 
much more information about the neutron.] A very useful index to study the neutron 
form factors is the electric square radius. The neutron electric square radius quantifies 
the slope of the Geu form factor at Q 2 = 0. 

In the present model we can quantify the contribution of the valence quarks and 
the pion cloud using the decomposition (|6.4[) and (|6.5I) for the neutron electric radius 
(X = E). In the present case the result is 

r| n = (-0.097) 6 + (-0.016) w fm 2 . (6.6) 

From this we conclude that the pion cloud gives about 14% of the total result 
[—0.113 fm 2 , very close to the experimental result —0.116 fm 2 ]. Again the pion cloud 
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0.643 
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0.001 
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0.325 


-0.005 
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0.382 


0.021 


0.403 


0.218 


0.050 


0.218 



Table 8. Octet baryon square radii in vacuum, and the bare and pion cloud 
contributions. 

contribution is still much smaller than suggested by an explicit calculation in the 
CBM [69]. Nevertheless, this is the best we have been able within the constraints of 
the fit to the whole octet. We note that the sign of the bare contribution is consistent 
with the lattice QCD data as we will explain next. Since Ge-u is positive and increases 
with Q 2 (positive derivative) it gives a negative contribution for r 2 En , according to the 
definition (|6.2p [with the replacement Geu(0) —> !]■ 

The effects of the valence quarks and pion cloud effects for the octet baryon 
radii were also estimated in (48] . There, the formalism used is the heavy baryon chiral 
perturbation theory applied to the lattice QCD data from [49]. In that work of Wang et 
al. [48] the quenched lattice results are corrected by subtracting the quenched and finite 
volume effects to extract the full QCD result in the physical limit of infinite volume. 
In this case we estimate the bare contribution from the quenched QCD contribution 
(QQCD), since in that work it is the component that better approximates the valence 
quark content. For simplicity we will ignore the uncertainties of the estimates. In 
that case pion cloud is the dominant contribution (67%) and the slope of the core 
contribution is positive. The last result suggests that the contribution of the quark 
core to Geu is positive. This is qualitatively consistent with our results, but it is 
important to quantify the contribution of each term. 

In our model the valence quark contribution is determined by the fit to the 
lattice QCD data and also the nucleon physical data. In both cases the neutron 
data are included. As discussed in the previous work [27] the spectator quark model 
simulates the magnitude and sign of the lattice QCD data [55J for Geu used in the 
calibration. That is a consequence of the VMD parametrization where there is an 
asymmetry between the up (it) and down (d) quarks electromagnetic structure. This 
is represented by a distinct parametrization of the isoscalar /i + and isovector /i_ 
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Q 2 (GeV 2 ) Q 2 (GeV 2 ) 



Figure 2. Nucleon electric form factors calculated in lattice, the present 
model (solid line), and the model from 27 (dashed line). For the physical point 
m w = 138 MeV, only the contributions from the core are included (without the 
pion cloud). 

quark form factors. In particular we have /i+(Q 2 ) > /i-(Q 2 ), leading to the resuh{jj| 
Geu(Q 2 ) > for lattice and also for the bare core in the physical case. We note that 

|| Using the expression for the nucleon Ge given by equation (28) of 20 , which is equivalent to the 
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Figure 3. Nucleoli form factors fitted at the physical point compared with the 
experimental data. The bands represent the pion cloud contributions. 



although quenched QCD simulations preserve the isospin symmetry, the full QCD 
simulations and the nature violate (in small degree) the isospin symmetry. 

In figure[2]we present the results of the model for the nucleon electric form factors 
in lattice (Ge), compared with the lattice data from [88] used in the present fit for 
the cases m, = 354, 495 MeV. The calculation is made using the model extended to 
the lattice regime based on the VMD parametrization [see [27] for more details] . The 
contributions of the bare core for the physical case (m-n = 138 MeV) are also presented 
and compared with the data. In figure [5] we present also the results of the previous 
model [27], which show that the both models have similar results. 

As we can observe from lattice data (m^ = 354 MeV, and m, = 495 MeV) for 
the neutron, the model results are close to the lattice data, although there are some 
overestimates of the data. Note also that the neutron lattice data are positive. Since 



one used in this work without pion cloud, we have 

(/l+-/l-)-4T^(/2+-/2-) 



Geti — —B 



where B is the overlap integral (normalized to 1 at Q 2 = 0). As the Pauli form factors are suppressed 
in the small Q 2 region, the difference between /i+ and /i_ determines the sign of Geu near Q 2 = 0. 
Taking Ge p — Gs n we get a dependence in the isovector form factors (i = 1,2). 
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Figure 4. Proton electromagnetic form factors calculated for p = (vacuum) 
with a decomposition of the valence and pion cloud contributions (left panel), 
and the ratios to those of the p = (right panel) for p = 0.5 po (dashed line) and 
1.0 po (dash-doted line). 



in this lattice QCD regime the pion cloud effects are not important, we interpret the 
result as a consequence of the isospin asymmetry of the model for the valence quarks 
as already mentioned. The figure supports our motivation to apply the quark model to 
the lattice QCD regime, using a model with isospin breaking. One can argue however 
that since we consider a quark current parameterized by a VMD structure, some pion 
cloud effects can also be in part of the current, particularly in the term associated with 
the p-pole. [For the high Q 2 region we can claim that the quark structure is the only 
one that survives consistently with the data.] To analyze better this point we look 
again for the neutron electric radius. Since the Geu data in lattice show a positive 
increasing function of Q 2 in the low Q 2 region, this implies a negative contribution of 
the bare core to the charge radius [see equation (l6.2[) ]. The simple comparison of our 
results with those of Wang et al. [35] based on the relative contributions mentioned 
before (67% of pion cloud in r 2 En ), can lead a conclusion that about 53% of our VMD 
parametrization may be pion cloud (67% from Wang's result, and minus 14% in our 
estimate). Notice, however, that Wang et al. uses a different structure for the pion 
cloud. Contrary to a more common representation for the pion cloud effects in the 
nucleon system, as the one we adopt (see figured]) following [SU] [HI], the work of 
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Figure 5. Same as the caption of FIG. |4] but for the neutron. We note 
that, as discussed in the text, some of the pion cloud contribution in models 
such as CBM |69l [70] are included through VMD in the valence terms in this 
parametrization. 



Wang et al. [48] also includes a diagram with a pion double vertex that comes from 
full QCD. This diagram can also be a source of the part for the ~ 50% difference 
between our estimate and that of [48]. Future studies using lattice QCD simulation 
data to reduce the model dependence for the pion cloud contributions, combined with 
a precise estimate for the valence quark structure in the intermediate Q 2 region, can 
help to pin down the effective contributions of these degrees of freedom. 

It can be useful for the nucleon case to compare the lattice QCD data from 
[55] . and our results using the other lattice QCD simulations. There are for instance 
different data in [92], from the same group, than the data we have used in our 
calibration. Unfortunately the lattice QCD studies of the nucleon form factors are 
mainly performed for the nucleon isovector form factors like — Ge p — Geu and 
G\j = G]\ip — GMn, and results for Geu have not been published in general. Examples 
are the works of QCDSF collaboration [93] and the Cyprus group [94] [95]. The main 
reason to avoid determining the form factors for proton and neutron separately, is 
because the isoscalar form factors, G E = Ge p + Geu and Gfj = Gm p + Gm u , should 
include contributions from disconnected diagrams [93], which is a complex task with 
the present lattice QCD resources. The results extracted from [49] [93] [96] are too 
imprecise to draw any conclusion about the magnitude and the sign of Geu on lattice. 
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Nucleon isoscalar form factors were calculated recently by the LHPC collaboration |97] . 
and they show also a positive result for Geu for pion masses similar to the ones in 
[88] . but larger in magnitude. 

Although there are other lattice simulations of the nucleon form factors, we keep 
our preference for the data in [88], even if the data are affected by some systematic 
errors as mentioned already [27j . The main reason is the fact that it is the only work 
extended to the octet baryons. Therefore, the limitations in the neutron form factor 
results can be compensated by the inclusion of the data for E + '~ and S 0, ~. 

The results of the fit for the nucleon form factors are presented in figure [3] and 
compared with the selected data EEl IZ3 EEl ESI UTS IHH IHH H31 HH IHH IH3 

(see [27] for a detailed discussion about the database). In figure [3] the solid line gives 
the full result (valence plus pion cloud) and the bands the effect of the pion cloud. 
The results for the remaining octet baryon members will be presented in the next 
section and extended to the nuclear medium. Recall that except for the magnetic 
moments there is no data for the A, E and S systems. Results in figure [3] show the 
dominance of the quark core and that the pion cloud effect is restricted to the small 
Q 2 region. As for the other octet members the form factors are mainly determined by 
the lattice QCD data. The fit to the lattice data provides a good description (small 
X 2 per data point) for the E systems but are not so good for the S systems. As for 
the neutral particles, A, E° and 5°, our calibration provide only a crude estimate of 
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Figure 7. Same as the caption of FIG. |4] but for E + . The experimental 
magnetic moment is shown by the (red) dot. 



the core effect, since these systems are not constrained by lattice data (except for S°). 
Therefore, the separation between the core and pion cloud should be taken with care 
particularly for Ge (very small for neutral particles), where small pion cloud effects 
and residual quark core effects cannot be distinguished with precision. 

Other properties of the octet baryons in vacuum are listed in Tables [7] and [8] 
From Table [7] we can conclude that the present model can describe fairly well the 
octet square radii and the magnetic moments, except for S°. From Table [8] we note 
the dominance of the quark core effects on the octet baryon radii. The exception is 
the electric square radii of the charge neutral particles, as mentioned already. A note 
about the large pion cloud contribution for the S _ magnetic square radii, is in order. 
The large pion cloud contribution for r 2 M ^_ is essentially a consequence of the large 
contribution of the term in Bi (photon-pion coupling contribution for Fib)- in fact 
the same effect appears for E + , but the final value is reduced by the larger magnetic 
moment (factor 2.28), in the definition of r\ iB [see equation (|6.3|) ]. 

Before extending the model for the in-medium regime, we discuss the sensitivity of 
the model to the input data. As already mentioned, we restrict the range of variation of 
r^ s _ to 0.610±0.015 fm 2 , instead of using the experimental result 0.61±0.15 fm 2 [67 . 
The main effect of this constraint is the reduction of the pion cloud contribution for 
r^L-y- to a value « 0.2 fm 2 , similar to that for the proton. Also the contribution for the 
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Figure 8. Same as the caption of FIG. [3] but for S°. 
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E + electric square radius, (less than twice the value for E~) is similar in magnitude. 
We can then conclude that the parameters of the model are very sensitive to the values 
of r E11 - and rjj S+ , and the respective contributions of the pion cloud. Thus, precise 
measurements of r 2 E ^ and J"|, s+ should be very useful to improve the quality of the 
model, and also to check if the assumption r EJ ,- ~ r Ep , is justified or not. 

A few notes about the quality of the results for the octet baryon results are in 
order. The quality of the fit is better for the systems with a strange quark, A, E, than 
the ones with two strange quarks, S. We recall that the kaon cloud contribution 
neglected for the S-baryon in the present approach, may be important since the 
contribution of the pion cloud is small. Therefore the prediction for £ are expected to 
be more reliable than that for H. Also, as no lattice data for A and E° are used in the 
calibration of the model, the results for the charge neutral particles, A, £°,£ , have 
to be taken with caution, particularly for the separation of valence and pion cloud 
contributions as discussed about the radii. Nevertheless, we present the results for the 
all charge neutral particles for completeness. Finally, we note that the predictions for 
the high Q 2 region has also to be taken with care, since the lattice data used in the 
calibration are restricted to Q 2 < 1.5 GeV 2 . The nucleon case is an exception (see 
figure [3]), since (physical) data are available for high Q 2 . 
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6.2. Octet baryon electromagnetic form factors in medium 

We consider now the octet baryon electromagnetic form factors in the nuclear medium. 
The formalism and the parameters necessary have already been presented in section 
[5] We calculate G* E and G* M for nuclear matter densities p — 0.5 po and p = po with 
p = 0.15 fm -3 . The modification of the form factors in the nuclear medium may be 
characterized by the in-medium modified masses given in Table [3] and the modified 
coefficients f3 B given in Table [4] 

The results are presented in figures I41TTT1 On the left panel in each figure the 
results in vacuum (p = 0) are presented, where also the valence quark contribution 
(dashed line) and the pion cloud contribution (dash-dot line) are shown. On the right 
panel, in-medium to vacuum ratios of the form factors are shown. For the charge 
neutral particles (n, A, E°, 5°), instead, the absolute values of G* E (Q 2 ) are shown, 
since G_e(0) = G* E (0) — and because the values are in general small for finite Q 2 
compared to G* M (Q 2 ). For these cases the vacuum values are presented with the thick- 
solid lines. The experimental magnetic moments in vacuum |67j are also shown for 
the cases B = A, S 0, ~ with the filled circles. 

From the figures it is clear that generally the valence quark contributions are 
more than 80% of the total contribution of each octet baryon form factor, and this 
is in agreement with the results of [37]. The exceptions are the electric form factors 
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Figure 10. Same as the caption of FIG. [4] but for S°. The experimental 
magnetic moment is shown by the (red) dot. 



of the charge neutral particles n, A, E° and H°. To estimate the pion cloud effects 
in the nuclear medium, we note that the ratio of the pion coupling constant in- 
medium to vacuum, g* BB i / 'g^BB' , given in Table [4] is at most 10%. As the factors 
/3b depend quadratically from the ratio g* BB > / g-^BB' , the effect of the variation in 
medium can be up to 20%. Since the largest case of the pion cloud contribution in 
vacuum electromagnetic form factors is about 20% (aside from the electric form factors 
of the neutral baryons), the total modification in medium associated with the pion 
cloud is at most 4% (20% of 20%). Thus, the pion cloud effects in the nuclear medium 
are essentially the same as those in vacuum, which are shown in the left panels in 
figures [4"l ll II Thus, the significant modification of the form factors in medium which is 
shown in the right panels in figures l4"l lll[ is according to the in-medium modification 
of the valence quark contributions, which are much more sensitive to the in-medium 
modification than those of the pion cloud. The strong sensitivity of the valence quark 
contributions is a consequence of the medium modification of the vector meson masses 
(modification of the VMD based quark current) and the baryon masses, which are in 
the radial wavefunctions. 

However, the smallness of in-medium change of the pion cloud contributions for 
the electromagnetic form factors (at most 4%), does not apply for the electric form 
factors of n,A, S° and 5°, since a priori these electric form factors are very small 
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G* E ps Q 2 near Q 2 = 0. In these cases it would be much informative to look the 
direct contribution of the pion cloud contribution in medium, not the ratio, since it 
can be dominant. In this way, we show the pion cloud contributions in vacuum and 
in medium in figure 1121 By comparing the scales we can conclude that the pion cloud 
contributions are small for the neutron and H , but they are the leading contributions 
for A and S°. For the neutron the contributions are about 0.01, to be compared with 
the total magnitude 0.06 (see figure[5]). As for the S°, the pion cloud contributions are 
negligible (compare the scales of G* E in fjgureUDland figure H2]l. which is a consequence 
of the small coupling constant (/?= ~ 0.04). Finally for A and S°, as shown in figures 
[B] and [51 the pion cloud gives dominant contributions in the region Q 2 < 1 GeV 2 , 
although the valence quarks become dominant for higher Q 2 . Thus, we conclude 
that, although the pion cloud gives dominant contributions for A and E°, the medium 
modification of the form factors according to the pion cloud is small. 

We look now for the global result in medium. For the electric form factors 
we notice that the ratio G* e /Ge decreases from unity for increasing Q 2 (except the 
discussions for the charge neutral particles) , which means that the electric form factors 
in medium decrease faster than those in vacuum. The medium effect is very small for 
S~ (ratio is almost constant) reflecting the lower sensitivity of the strange quark to 
the medium modifications). As for the magnetic form factors they are enhanced by 
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the medium for small Q 2 . The effect is larger for larger densities. In [10] proton 
electromagnetic form factors in medium were studied. The results presented in jlOj 
show similar trends to those observed in the present study. Namely, the in-medium 
magnetic form factor is enhanced, while that of the electric form factor is quenched, 
where these results can explain the modification of the bound proton form factors 
measured at Jefferson Lab [5] . The modifications of the nucleon electromagnetic form 
factors in medium will be discussed in next subsection. 

In the literature we find only |98j studied medium modification of the octet baryon 
magnetic moments aside from the nucleon, using a different kind of quark-meson 
coupling model from the present one, "QMC" and modified quark-meson coupling 
model (MQMC). In their treatment, the pion or meson cloud is not included in 
the electromagnetic currents. (We discuss the nucleon case in next section.) They 
compared with "QMC" and MQMC results for the magnetic moments in nuclear 
medium. In their "QMC" , the in-medium magnetic moments for the A and S _ 
decrease compared to those in vacuum, although other octet magnetic moments are 
all increased in medium. This feature is different from the present approach, where all 
the octet baryon magnetic moments {G* M {Q 2 ) at Q 2 = 0) are enhanced in medium. 
In their MQMC, on the other hand, all the octet magnetic moments are enhanced, 
and this feature is the same as that of the present approach. 
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6.2.1. Nucleon electromagnetic form factors in-medium. We study now in detail the 
medium modification of the nucleon electromagnetic form factors. A very interesting 
quantity is the ratio of the electric to the magnetic form factors, 

Rn = £K (6-7) 

"J/JV 

which can be also calculated in vacuum (denoted by Rn)- The study of this ratio for 
the proton in vacuum serves as a fundamental quantity in the present understanding of 
the proton structure, as was measured at Jefferson Lab, that R p significantly deviates 
from a constant [74] E5] • Similar studies were also made for the neutron [81] [82] [99] . 

In figure [13] we present our predictions for the both, proton and neutron ratios 
in medium, for nuclear densities p = 0.5 po and p = 1.0 po- Note that p = (vacuum) 
case was already compared with the experimental data in the present model since 
it fits both Ge and Gm (see figure [5]). Experimentally, direct access for the ratio 
R* N in-medium seems to be impossible at present. However, we can get indirect 
information for the in-medium ratios using the results of proton (nucleon) recoil 
polarization experiments in nuclei by the measurement of the polarization-transfer 
super-ratio for the nucleon, 

K=p^M. (6.8) 

The experiments were performed for the proton using 1 H, 4 He, and 16 targets [HG3IH], 
and also planed for the neutron case [9] . 

The results of our predictions for 1Z for the proton and the neutron arc 
presented in figure 14. Some other calculations for the proton case can be found 
in [TU] ED E3 EM EM EM EM EU, and summarized in For a detailed 

discussions about the super-ratio of the proton and the neutron, see M [5] 1105] . 

First, we discuss the proton super-ratio shown in the left panel in figure [141 with 
the data with the 4 He target from [4] [5] HjJ. Our results for p = 0.5 po (po = 0.15 
fm~ 3 ) reproduce the data. While the average nuclear density of 4 He is expected to be 
slightly higher (0.74 po in a QMC estimate) the effect of absorption means that the 
active nucleon tents to be in the nuclear surface. Thus, the present model predicts 
the observed trend of the reduction for the super-ratio. 
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Next, turning the discussion to the neutron case, our results predict an 
enhancement of the super-ratio, contrary to the case of the proton. The enhancement 
was also predicted by Nambu-Jona-Lasinio (NJL) model and relativistic light front 
constituent quark model (LFCQ) |S]. However, the difference is the Q 2 dependence. 
In the case of the NJL model, the super-ratio monotonically decreases with increasing 
Q 2 , while it stays almost constant in LFCQ. The present result shows appreciable Q 2 
dependence, namely the ratio increases up to around Q 2 = 0.3 GeV 2 , and gradually 
decreases with increasing Q 2 . Thus, our model predicts that the modification can 
be maximally observed around Q 2 = 0.3 GeV 2 . This point may be taken into 
consideration in the experiments planned [9]. 

7. Summary and conclusions 

In this work we have presented a model for the octet baryon electromagnetic form 
factors in vacuum and in the nuclear medium. The model is based on a constituent 
quark formalism but includes also a phenomenological parametrization for the pion 
cloud motivated by xPT. The octet baryon and relevant meson properties in the 
nuclear medium (masses and coupling constants) are determined using the QMC 
model. The effects of final state interactions and meson exchange current are not 
explicitly included. 

The model is calibrated by the octet baryon lattice data (electromagnetic form 
factors), as well as physical data, such as the nucleon electromagnetic form factors, 
octet baryon magnetic moments and the available octet baryon radii data. Lattice data 
give stringent constraints for the valence quark structure of the octet baryons. The 
remaining data combined with the estimate of the valence quark core contributions, 
constrain the parametrization of the pion cloud. 

The fit is very sensitive to the neutron lattice data for the electric form factor and 
also to the value of the S _ electric radius. The neutron data, lattice and physical, 
are well described by a valence quark model with isospin breaking. An accurate 
description of the data, including the octet radii (nucleon and E _ ) is also obtained 
with a small value for electric square radii for E~ (w 0.6 fm 2 ) which is a consequence 
of a small pion cloud contribution for the electric radius (same order as that for the 
proton). 

Future improvements are possible. Accurate lattice data, particularly for 
the neutron electric form factor (Geu) can clarify the role of the valence quark 
contributions, and eventually demand a refit of the quark current used in this work. 
The available lattice data support a quark electromagnetic current with an isospin 
breaking, but new data with a smaller effect for Geu may require a model with 
almost no isospin breaking. The model can also be improved by including a pion 
cloud parametrization derived from first principle QCD, utilizing the results of lattice 
simulations (see Ref. [IE])- 

The explicit inclusion of meson exchange current corrections in the electromag- 
netic form factors in a consistent manner, may be appreciable at high momentum 
transfer, especially if the absolute values of form factors became very small, or for the 
electric form factors of neutral baryons. 

We predict that all the octet baryon magnetic form factors in medium will be 
enhanced in the small Q 2 region, and the magnetic moments of the octet baryons in 
medium are enhanced. The model predicts that Q 2 dependence of the octet baryon 
electric form factors in medium decreases faster than those in vacuum. Furthermore, 
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Figure 14. Super-ratios in nuclear matter calculated for the proton (left) and 
neutron (right). Data for proton are from 4 5] [6]. 



the present model predicts also a decrease of the super-ratio for the proton in the 
nuclear medium, which is in agreement with the observed results. On the other hand, 
for the neutron super-ratio in the nuclear medium, the present model predicts the 
enhancement, which has its maximum around Q 2 — 0.3 GeV 2 , and this may be useful 
information for the planned experiments. 
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